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Tackling the infrared problem of thermal QCD mit-ctp-3421 
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Perturbative calculations of corrections to the behavior of an ideal gas of quarks and gluons, the limit that 
is formally realized at infinite temperature, are obstructed by severe infrared divergences. The limits to com- 
putability that the infrared problem poses can be overcome in the framework of dimensionally reduced effective 
theories. Here, we give details on the evaluation of the highest perturbative coefficient needed for this setup, in 
the continuum. 
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1. INTRODUCTION 

The theory of strong interactions, Quantum 
Chromodynamics (QCD), is guaranteed to be ac- 
cessible to perturbative methods once one of its 
parameters, the temperature T, is increased to- 
wards asymptotically high values. This general 
statement relies solely on the well-known prop- 
erty of asymptotic freedom. 

In practice, however, calculations of corrections 
to the behavior of an ideal gas of quarks and glu- 
ons, the limit that is formally realized at infinite 
T, are obstructed by severe infrared (IR) diver- 
gences for every observable, there exists an 
order of the perturbative expansion to which an 
infinite number of Feynman diagrams contribute. 
No method is known how to re-sum these infinite 
classes of diagrams, a fact that seriously obstructs 
progress in the field of thermal QCD. 

It is known how to evade this obstruction us- 
ing dimensionally reduced effective theories. The 
key idea is to map the infrared sector of thermal 
QCD onto a three-dimensional pure gauge theory 
1 112131 1J . whose contribution, being a pure num- 
ber, could be extracted numerically by Monte- 
Carlo simulations. While the expansion of the 
QCD pressure in the effective theory framework, 
up to the order where IR contributions are rele- 
vant, is now known analytically 0], realizing the 
numerical extraction of the yet-unknown number 
emerging from the IR sector is a challenging open 
problem, with the main complication that high- 
order matching between lattice and continuum 
regularization schemes is necessary [S]- 



2. SETUP 

Let us now switch gears and focus on one of 
the main building blocks of the procedure, while 
for a detailed description of the setup as well as 
notation and further references, we refer to 0]. 
In particular, we want to compute the (negative) 
3d vacuum energy density of a pure SU(N) gauge 
theory, 
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which in a weak-coupling expansion can be writ- 
ten as the sum of all connected vacuum graphs 
containing gluons and ghosts. Since the theory 
is confining, the computation involves IR diver- 
gent integrals (starting at the 4- loop level here), 
forbidding a perturbative evaluation of the full 
vacuum energy. One can however obtain its log- 
arithmic ultraviolet divergence. 

Note that in 3d the coupling constant g is di- 
mensionful, hence the full answer must be of the 
form 

a G (- + 8\n-^—)+j3 G +0(e)] , 
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CaQ 1 is a dynamically generated in- 
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frarcd scale in the confining theory, and Ca = N 
and (1a = N 2 — 1 are the Casimir and the dimen- 
sion of the adjoint representation, respectively. 
Because of super-rcnormalizability, the coefficient 
og can then be computed in 4-loop perturbation 
theory, even if the constant part (3q cannot. 

If we just carry out the 4-loop computation in 
strict dimensional regularisation, the result van- 
ishes, because there are no perturbative mass 
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scales in the problem. This means that UV and 
IR divergences (erroneously) cancel against each 
other. Therefore, we have to be more careful in 
order to determine cue ■ To regulate the IR diver- 
gences, we introduce by hand a mass scale, m 2 , 
into the gauge field (and ghost) propagators. One 
has to keep in mind, however, that now only the 
coefficient ac multiplying 1/e is physically mean- 
ingful, as it contains the desired gauge indepen- 
dent ultraviolet divergence. On the contrary, the 
constant part depends on the gauge parameter 
£, because the introduction of m 2 breaks gauge 
invariance, and has nothing to do with [3q. 

Note that e.g. diagrams with self-energy inser- 
tions can have IR sub-divergences, since IR diver- 
gences are known to be present in the 3d 2-loop 
gluon propagator. To avoid the problem of over- 
lapping IR divergences from the outset, we have 
hence chosen to employ the mass parameter rigor- 
ously, i.e. by rewriting every 1/p 2 as l/(p 2 + m 2 ). 

This leaves us within the class of fully mas- 
sive integrals. The computation can be divided 
in three parts. Roughly, those are (1) diagram 
generation |H] , specification of Feynman rules and 
color algebra, (2) reduction to master integrals 
|7I8| . (3) expansion in d = 3 — 2e dimensions. 

We will refrain from commenting on the first 
two parts of the computation here, since they are 
well documented in the references given above. 
Due to the complexity of the computation, both 
steps are automatized, allowing for the handling 
of a large set of diagrams. 

3. MASTER INTEGRAL REPRESEN- 
TATION 

Let us now give a little more detail on part (3) 
of the computation. At this point, all diagrams 
are expressed in terms of 19 scalar master inte- 
grals, which are enumerated in [Hj. The general 
structure is 

where d is still an arbitrary (space-time) dimen- 
sion. Only now do we need to specify d = 3 — 2e. 

While it is trivial to expand the polynomial 
prefactor in e, considerable effort has to be put 



into obtaining the expansion for the master inte- 
grals to the depth required. Since we need the 
e-poles only, it would seem sufficient to compute 
the divergent parts of all master integrals. It 
turns out, however, that the prefactor develops 
poles as well around 3 dimensions, having terms 
proportional to 1 / (d — 3) multiplying 10 of the 
master integrals, and even double poles in 4 of 
those cases. 

A crucial simplification can be made by exploit- 
ing the freedom of choosing the basis of master in- 
tegrals to represent the sum of diagrams Eq. J5J. 
Going back to the tabulated relations between in- 
tegrals that were derived by partial integration 
and used in part (2), we found two most useful 
relations: 

(3D - -^02) 
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(2d- 7)(2d- 5) 
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Notation: each line represents a massive scalar 
propagator, a dot on a line means an extra power, 
vertices have no structure. Trading the two mas- 
ter integrals on the lhs of the above equations 
for the first ones on the rhs respectively (all oth- 
ers are already included in the basis), the d- 
dimcnsional representation Eq. of course still 
holds, albeit with a 'primed' version of the basis, 

^jf^^^^M^d). (5) 
(4^)4 z-. poly-(d) 

In this new basis, none of the prefactors has a 
double pole in 3d, while only 7 members of the 
new 'primed' basis are multiplied by a single pole. 
It is not excluded that there exists a choice of 
basis for which the prefactors never get singular, 
but this choice is currently not known to us. 
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4. EXPANSION 

It turns out that (almost) all integrals are 
known analytically to the order needed for obtain- 
ing the poles in the sum of all diagrams. Lower 
loop cases have been treated in ^U|, while ana- 
lytic results for the divergences of all 3d 4-loop 
master integrals as well as numerical and some 
analytic results for their constant parts as well as 
the 0(e) term of the 2-loop sunset integral can 
be found in By an amusing relation specific 
to 3d, namely the fact that the leading term of 
the 3d 1-loop scalar 2-point integral is an arctan, 
whose derivative with respect to a mass looks like 
a propagator with double mass, it is furthermore 
possible to relate the leading term of one of the 
4-loop master integrals to a 3-loop case |12| : 



(6) 



There are however 2 master integrals (out of 
the 7 which get multiplied by a 1/e from the pref- 
actor) whose constant term we do not yet know 
analytically. Let us denote their leading parts by 
Xi and £3 (by naive power-counting, it is easy to 
see that both are UV finite), 
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2:3 + 0(e) (7) 



Filling in the known expansions for the mas- 
ter integrals as well as expanding the prefactors, 
higher poles cancel in the sum of diagrams, and 
we are left with a single pole only: 
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The polynomial p is of order 6 in the gauge pa- 
rameter £ and contains, besides a collection of 
numbers like 7r 2 , In 2 and dilogarithms, the two 
unknowns x-i and X3. Clearly, in order for the 
result to be gauge independent, all ^-dependence 
has to vanish once x-i and X3 are known. We can 
now reverse the argument and try to fix these 
constants by requiring gauge independence. In- 
specting the polynomial, it turns out to have a 
very simple structure: 



p(0 = a G + (x2 - 6x3 - b) a C 
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a G 
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90 
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6144 



0.195715.. 
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where the Cj are pure numbers and b = \A 2 \ 
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2(ln2) 2 - |(ln3) 2 - (ln5) 2 - 
We have checked by 



21n21n5 + 31n31n5 + \ . 
numerical integration that 

£2-6x3 = b « -0.00200966335 .. . (11) 

to nine significant digits, hence establishing 
Eq. (| 1 U|) as our main result for the logarithmic 
divergence of 3d pure gauge theory. 
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